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AB initio quantum chemical techniques have been used to 
investigate weakly bound complexes of H 2 0 and S0 2 . An energy 
gradient program was used to locate stable structures for the 
H 2 0 J S0 2 complexes and SCF calculations were carried out to 
determine the binding energies of complexes with multiple water 
molecules. A 4-31G basis set was used for most potential energy 
searches. More accurate basis sets including a generally contracted 
basis set with d orbitals on the sulfur were used for geometry 
and binding energy verification. For single water complexes, 
five different stable geometries were located with binding energies 
between 4 and 11 Kcal mol" 1 ^ suggesting a binding shell for H 2 0 
around S0 2 and a mechanism for the formation of an S0 2 -containing 
water droplet. Calculations on one of the complexes utilizing a 
larger double zeta basis and d functions on the sulfur atom lead 
to adjusted binding energies in the range 3 to 8 Kcal mol -1 . 
Very little charge transfer between S0 2 and H 2 0 was present. 
Addition of more than one H 2 0 was found to be energetically 
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favorable although the addition of the fourth water in certain ; 

— / 

geometries did not increase the stability of the complex. An [ 
alternative mechanism for the troposheric, gas phase production ; 
of acid rain is suggested. 

Introduction 

The oxidation of SO 2 is an important chemical process. The 
phenomenon of "acid rain" involves the oxidation of sulfur 
containing species including SO 2 to form SO 3 . Hydration is 
involved to give H 2 S0 4 which in turn is ultimately precipitated 
as sulfates . 1 The influence which water plays on the kinetics of 
oxidation is less clear. However, Reiss noted a 4.7 x 10 3 fold 
increase in the rate of photo oxidation of S0 2 in the 283% 
relatively humidity of a Wilson cloud chamber over that in the 
absence of water 2 and Ravishankara noted that the formation of 
aerosol upon the irradiation of gaseous mixtures containing H 2 0 
and SO 2 and the formation of aerosol in H 2 O • S0 2 * O 2 mixtures. 
Ravishankara further noted that aerosol formation occured only 
for sufficiently high water vapor concentrations . 3 

The stability of a complex of SO 2 and H 2 O was first suggested 
by Phillips on the basis of results of semiempirical calculations . 4 
Holland and Castleman using CNDO/2 calculations suggested that 
an SO 2 * H 2 O adduct should be stable in the gas phase and reported 
a binding energy of 145 mH and a large dipole moment.^ 



Evidence for the existence of H 2 0 • S0 2 produced from the 
photolysis of H 2 S in solid 0 2 at 15K by ultraviolet light is 
presented by Tso and Lee who conducted an FTIR study of the 
photolysis mixture . 6 Identification was based on shifted vibrational 
bands of S0 2 and H 2 0. 

The stability of S0 2 -water complexes has also been suggested 
from molecular beam measurements, and electric deflection experiments 
indicate a large dipole moment for the complex. 

Ab initio calculation on S0 2 • NH 3 8 and S0 2 • HF 9 show stable 
complexes for a variety of geometries. An important pair of 
calculations was done on the closely related complex H 2 0 • S0 3 and 
its rearrangement to H 2 S0 4 . A semiempirical CNDO/2 calculation 
by Holland and Castleman -10 has been followed by a recent £b initio 
calculation by Chen and Plummer . 11 After the completion of the 
calculations reported here, the authors became aware of a articles 
by Plummer 1 ^ and by Plummer, Chen, and Law 11 which reported ab initio 
results on two structures of the monohydrate of S0 2 obtained using 
an ST0-3G basis set augmented by a set of d functions on the 
sulfur atom. Plummer reported binding energies of 10 and 2.3 
Kcal/mol respectively for the two monohydrate structures studied. 
The structures reported by Plummer are similar to two of the 
monohydrates investigated in the present study , but the results 
differ significantly in detail from those reported here (vide infra) . 
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In an effort to shed some light on a possible role of H 2 0 in 
the oxidation of S0 2 and elucidate the structure and energetics 
of formation of water complexes of S0 2 , we have undertaken an a ]2 
initio study. 

S0 2 • (H 2 0) ! 

To explore the possibility of bound conformations of H 2 0 • 
S0 2 an SCF energy gradient program 14 was used to locate stable 
structures for the complex. The 4-31G basis set was used for 
the initial potential energy searches. The minimum energy 
geometry was established when the largest component of the energy 
gradient was reduced to the order of 1 milliHartree per Bohr. 

Five distinct conformations of the single water complex were 
found to be bound. These are shown in Figure 1. The total 
energies for the complexes and the molecular fragments as well 
as the binding energies are given in Table 1. The minimum energy 
geometry information is given in Table 2. 

Conformations III and IV bear the closest relationship to 
classical hydrogen bonded structures, having respectively, one 
and two such bonds, and no significant binding interaction 
between the water oxygen and the sulfur atom. Conformation I, 
the most weakly bound, is at the opposite end of the scale with 
no hydrogen bonding, but has a favorable dipole interaction and a 
loose coordination between the water oxygen and the sulfur atom. 
Conformations II and V exhibit intermediate characteristics with 
both hydrogen bonding and 0-S coordination. Except for the 
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interplanar angle between the two moieties, conformation V, the 
most strongly bound, exhibits the most favorable alignment of the 
fragment dipoles. 

Hydrogen transfer with subsequent charge and structural 
relaxation from either conformation II or V would lead to the 
presumed structure of sulfurous acid. The coordination of the 
water oxygen and sulfur atom in conformations I, II, and V could 
also facilitate other reactions leading to oxidation of SO 2 • 

In conformation III the out-of-plane angle of rotation of 
the nonhydrogen-bonded OH is very soft energetically and no 
angle is favored by more than a few tenths of a milliHartree. 
Conformation III is similar to the more weakly bound structure 
reported by Plummer except that the present conformation is for a 
cis arrangement of the two fragments while Plummer's structure is 
trans. The interfragment separations resulting from the two 
investigations were almost the same but Plummer found a binding 
energy (2.3 Kcal/mole) which is significantly smaller than that 
given in Table I for conformation III. The present structure 
provides a favorable electrostatic interaction between H2 and 01 
of conformation III (Fig 1) and would be expected to yield a 
lower energy. Although the potential energy was found to be 
weakly dependent on the rotation angle of the nonhydrogen-bonded 
OH near the cis minimum, the calculations along this coordinate 
were not carried out at sufficiently large angles to determine if 
there is a significant potential energy barrier between the cis 
and trans structures. 
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Conformation V was suggested by the results of a semiempirical 
calculation 5 but the minimum energy geometry result here indicates 
a preference for a cis arrangement rather than the trans arrangement 
of the semiempirical calculation. The potential energy surface 
search which led to conformation V was started near the geometry 
found by Castleman and migrated toward the cis conformation as 
the iterations continued. However, the calculations were converging 
quite slowly and were terminated with a structure similar to the 
more stable of the structures reported by Plummer. The calculations 
were restarted in a cis configuration and subsequent iterations 
produced a fairly rapid decrease in the interfragment bond length 
and a significant increase in the binding energy. In light of 
Plummer's results it appears that a trans structure was probably 
skipped over in carrying out the search. Although the literature 
contains a significant amount of information on the performance 
of various basis sets on covalently bonded molecules, the information 
available for evaluating the results obtained on weakly bound 
complexes using different basis sets is much more limited and 
somewhat contradictory. For example, while Raffenetti and Phillips 15 
found that SCF calculations utilizing an unpolarized ST 0 - 3 G basis 
set overestimated the binding energy of H3N-HCI by almost a 
factor of two, Lucchese, Haber, and Schaefer 8 found that this 
basis set underestimated the binding energy of H3N-SO2 by a similar 
factor. The results obtained by the latter authors and the 
results of the present investigation are in agreement with respect 
to the relative importance of the completeness of the s and p 
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basis and the inclusion of d functions on the sulfur atom for the 
determination of binding energies. The fact that Plummer, Chen, 
and Law 13 found that the binding energy of the structure reported 
by Plummer 12 decreased by 0.8 Kcal/mole when the flexibility of the 
STO-3G* basis was increased somewhat by splitting it to 3-21G* 
suggests that the minimum basis set results in an overestimation 
of the binding energy in this case. Qualitatively, the increased 
interaction energy resulting from the favorable arrangement of 
the fragment dipole moments in the cis structure also suggests 
that this structure should be more stable than the trans structure. 
The shorter interfragment bond length found for the cis conformation 
would be expected for a more tightly bound structure. 

Figures 2 in the articles by Plummer and Plummer, Chen, and 
Law give a value of 2.7 Debye for the dipole moment of their more 
stable monohydrate structure 12 ' 13 . The value obtained for the cis 
structure in the present treatment with the DZP basis set is 
slightly smaller than these values, but the trans structure would 
be expected to have a larger dipole moment them the cis conformation. 

To examine the effect of basis set size on the binding 
energy and geometry two additional SCF energy gradient calculations 
were performed on conformation V. In the first calculation a 
double zeta basis set was employed. The primitive basis set, a 
Huzinaga sulfur (12s/8p) , oxygen (9s/5p) , hydrogen (4s) set, was 
segmentally contracted. 16 ' 17 A Dunning contraction of the hydrogen 
and oxygen orbitals, 0 (9s/5p)/<6111/41>, H(4s)/<31>, and a McLean 
contraction of the sulfur orbitals, S (12s/8p) /<62211/621> , were 
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chosen . 18 ' 19 The binding energy decreased by 3.0 mH, a rough measure 
of the superposition error in the poorer basis set. In the 
second calculation a set of sulfur 3d orbitals ( 33 = 0.6) 
was added to this basis set. The polarization functions reduced 
the S0 2 bond lengths and bond angle, decreased the S 0 2 dipole 
moment, and lowered the binding energy of the complex by an 
additional 2.4 mH. The results of these two calculations are 
given in Tables 1 and 2. 

Although the changes in the S0 2 structure and the complex 
binding energy were significant, the changes in the interfragment 
structural parameters were small. The decrease in the binding 
energy in the going from the 4-3 1G basis set to the more complete 
basis is very close to the percentage reduction which has been 
observed in previous treatments of similar complexes . 15 The 
present results indicate that more than half of the discrepancy 
in the binding energy is due to the incompleteness of the s and p 
functions in the poorer basis with the balance being due to the 
absence of polarization functions. Further augmentation of the 
basis set with polarization functions would be expected to lead 
to monomer structures and dipole moments in better agreement with 
experiment and a slight further decrease in the predicted binding 
energies. The inclusion of electron correlation corrects for the 
tendency of the SCF model to predict too much charge separation 
and produces more accurate dipole moments for monomers. The net 
effect of including correlation for complexes such as those 
investigated here is to increase the binding energy of the complex 
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by an increment on the order of 1 Kcal/mol . 15 > 20 < 21 
S0 2 • (H 2 0) n 

Selected two-water complexes were studied using the energy 
gradient program with the 4-3 1G basis set. The conformations 
are shown in Figure 2 . Conformations VI and VII are mirror 
image doublings of conformations I and V, respectively. 

The interatomic distances did not change significantly in 
going to the two-water complexes. The energies are given in 
Table 1. In both cases the binding energy of the two-water 
complexes are a little less than twice the binding energies of 
the one-water complexes indicating relatively little interaction 
between the water molecules and almost independent addition of 
the second water molecule. 

In order to examine the energetics of adding multiple waters 
to S0 2 with the double zeta plus sulfur d orbital basis set, SCF 
calculations were run using a generally contracted basis set at 
fixed geometries predicted by the earlier gradient results. The 
primitive basis set was again the Huzinaga sulfur (12s/8p), 
oxygen (9s/5p) , hydrogen (4s) set generally contracted to <4s/3p>, 
<3s/2p> , and <2s>, respectively. 16 ' 17 ' 22 This set was augmented 
by a single set of 3d functions ( 3d = 0.6) on the sulfur atom. 
This basis set was similar to the previous double zeta plus d 
set and the binding energy of conformation V obtained with this 
new basis set differed by only 0.2 mH from the previous result. 
Based on the results of optimizing the structure of complex V, 
the major effects of relaxing the geometry rather than using the 
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fixed geometries would be expected to be changes in the S0 2 
structure with some attendent lowering of the energy of the 
complex. Thus the binding energies obtained using fixed geometries 
are lower limit SCF binding energies for the complexes. 

The multiple-water complexes studied are shown in Figure 2. 
The total energies and binding energies are given in Table 1. 
Conformation I with this basis set produces a binding energy 
which is 72% of the binding energy derived from the earlier 
4 — 3 ig calculation. This reduction in binding energy is virtually 
identical to that obtained for conformation V. The multiple-water 
complexes S0 2 (H 2 0) n with n = 3 and 4 differ by only a few mH 
from the sum of the binding energies of the separate bound 
pairs. This assumes that the energies of the conformations II 
and IV (which were not calculated with the larger basis) can be 
approximated by taking 70% of the binding energy obtained with 
the 4-31G basis set. The binding energy of the five-water complex 
is, however, significantly less than the sum of the binding 
energies of the separate bound-pairs. 

The implications of the additivity of the binding energies 
supports the notion that the bonding interaction is nearly 
purely electrostatic. Only when the waters were close enough to 
interact with each other, as in the five-water complex, did the 
binding energy fall significantly below the sum of the binding 
energies of the separate bound pairs. An interesting corollary 
to the approximate additivity of the binding energies occurred 
in the three water complex, conformation VIII, which was slightly 
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more stable than predicted from the binding energy additivities 
perhaps due to a favorable alignment of the water dipoles. 
These results suggest that the binding energy of some multiple 
water complexes which could be modeled by superposition of smaller 
complexes, but which have not been investigated here, would 
deviate significantly from additivity due to water-water inter- 
actions. 

Although a Morakuma component analysis 23 of the energy components 
of the H 2 O - SO 2 interaction was not made, the results obtained 
here suggest that the analysis would follow that for the NH 3 • 
SO 2 complex 7 as discussed by Kollman. 2 ^ Mulliken population 
analyses 25 of the results of SCF calculations on the H 2 0-SO 2 
complex using the best basis set show a small charge transfer 
resulting in a small negative charge on the S0 2 . See Table 1. 
With the 4-31G basis set the charge transfer was about fifty 
percent larger, a consequence of basis set superposition error. 
The 4-3 1G basis set calculation for conformation IV produced the 
unique result of a charge transfer in the opposite direction 
from all the other results. Since the charge transfer and 
polarization affects are small, the minimum energy geometry is 
determined by a balance of the electrostatic energy and the 
exchange repulsion energy. However, unlike the NH 3 • S0 2 complex, 
the minimum energy orientation for the H 2 O and SO 2 fragments in 
conformation V is closer to that predicted solely from consideration 
of the dipole-dipole interaction. 
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The H-bonded complex, conformation III, has a binding energy 
(4.7 Kcal) which is virtually the same as that obtained for the 
HF • S0 2 complex (5Kcal) . 9 The minimum energy complex, conformation 
V, has a binding energy (10.9 Kcal) very close to that obtained 
for the NH 3 • S0 2 complex (9.3 Kcal). 8 It should be noted that 
inclusion of d orbitals in the sulfur basis set reduced the 
binding energy of NH 3 • S0 2 from 10.4 Kcal to 9.3 Kcal, an effect 
that was also observed in the H 2 0 • S0 2 calculations. The CNDO/2 
results of Castleman gave a much larger binding energy (98.3 
Kcal) and a different minimum energy geometry. 5 

In the analysis of his molecular beam electric deflection 
results Castleman, remarks on the existence of species with 
large dipole moments that he attributes to H 2 0 • S0 2 complexes. 7 
The dipole moments for some of the complexes were calculated in 
this study and are reported in Table 2. These values are all 
overestimated, particularly those from calculations which did 
not include sulfur d orbitals in the basis set. 

The gradient SCF program can also be used to produce force 
constants and fundamental frequencies for vibrational modes. 
However, the equilibrium geometries were not converged to sufficient 
accuracy to enable predicting the sizes and signs of the frequency 
shifts of S0 2 and H 2 0 vibrations upon complexation. Better 
convergence would allow comparison with the experimental IR 
spectroscopic results of Tso and Lee who claim to have evidence 
for the existence of the H 2 0 • S0 2 complex. 8 
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The gas phase reactions of S0 2 with a number of oxidants are 
exothermic, but the measured bimolecular rates for most such 
reactions are quite small. 26 In some cases, such as the reaction 
between S0 2 and CH 3 0 2 , complexes which subsequently decompose to 
the reactants are believed to be formed in preference to the 
bimolecular reaction. 26 Thus oxidants like 0 3 , N0 2 , H 2 0 2 , H0 2 , 
and CH 3 0 2 which have the highest concentrations in the trosposhere 
aren't considered important in the homogeneous oxidation of S0 2 , 
except in moderately to highly polluted air. Instead, most of 
the gas phase oxidation of S0 2 to H 2 S0 4 is believed to proceed by 
the following mechanism: 


03 
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S0 3 + H 2 0 + M— ?'h 2 S0 4 + M* (5) 

At 300K, this mechanism requires a steady state OH concentration 
of about 10 7 cm -3 in order to yield an SO 2 lifetime on the order 
of hours. Unfortunately, the distribution of OH in the troposphere 
is not known accurately. Since the above mechanism is driven by 
photochemistry and the lifetime of the OH radical is very short, 
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it is operative only in daylight hours. 

Heterogeneous removal of S0 2 by H 2 0 2 , N0 2 , and other species 
has been observed in aqueous solutions, however 1 ' 27 . Hetero- 
geneous oxidation of S0 2 in cloud droplets is now considered to 
be significant but remains to be accurately quantified. Estimates 
of S0 2 removal rates due to heterogeneous oxidation by H 2 0 2 
dissolved in cloud droplets range up to several hundred per cent 
per hour. 1 

Penzhorn and Canosa 28 have reported water vapor catalysis of 
the reaction with N0 2 . The results of the present investigation 
raise the possibility that complexes of S0 2 and H 2 0 may provide a 
homogeneous mechanism for catalyzing the oxidation of the former 
by oxidants which are present in the atmosphere at concentrations 
significantly greater than those of OH. This mechanism would not 
be subject to the pH quenching observed for many of the heterogeneous 
reactions. Such a mechanism could be important if the product of 
the equilibrium constant for the complex, the concentration of 
water vapor, the reaction rate, and the concentration of the 
oxidant were on the order of 10” 8 . This is equivalent to the 
oxidation rate due to 10 8 OH molecules cm” 2 at 300K. A somewhat 
smaller product would be significant for oxidants which do not 
drop to zero concentration in the absence of sunlight. 

Determining an accurate first principals estimate of the 
temperature dependent equilibrium constant would require knowledge 
of the barriers separating the various structures identified for 
the single water complex, barriers to free internal rotation, 
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vibrational frequencies, and the vibrational anharmonicities. 

A rough estimate of the range into which the equilibrium constant 
might be expected to fall can be obtained by comparison with the 
information available on other weakly bound complexes. Curtiss, 
Frurip, and Blander have determined K for both H 2 0 and D 2 0 dimers 
between 358K and 387K. 29 A 358K their results for D 2 0 give K as 
7 x 10 -19 cm 3 molec -1 . Extrapolating their results to 300K 
yields an association constant of about 1.4 x 10 -18 cm 3 molec -1 . 
The calculated electronic binding energy for H 2 0-S0 2 is higher 
than that of the D 2 0 dimer which would lead to a somewhat higher 
enthalpy of formation. Curtiss et. al obtained a standard 
entropy change of association of -18.67 cal/deg-mol for the D 2 0 
dimer. The possibility that the multiple binding sites for H 2 0 
in the H 2 0-SO 2 complex may contribute to the entropy change makes 
estimates of the latter subject to considerable uncertainity , but 
for the present purposes we may assume that it is approximately 
equal to that of the D 2 0 dimer. This would lead to a lower limit 
of about 10 -18 cm -3 molec -1 for the association constant with an 
upper limit as much as a factor of 100 larger. 

For the smaller of these two limits and water vapor 
concentrations appropriate for humid summer conditions, the reaction 
rate required to meet the above criteria would be on the order of 
10 -17 cm 3 molec -1 s -1 for ozone and 10 -13 cm 3 molec -1 s -1 for oxidants 
present at relative concentrations of a few parts per trillion. 
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For a number of oxidants, such as H2O2, which are present in 
moderately polluted air at relative concentrations on the order 
of 1 ppb, an intermediate rate of about 10”^® cm 3 molec ^s ^ would 
be required. With an equilibrium constant near the upper limit, 
the rates required for oxidants to be effective according to the 
criteria above would be proportionately lower. Currently recommended 
upper limits for the uncatalyzed reactions of trace oxidants 
other than OH with S 0 2 range from 5 x 10 -17 cm 3 molec” 1 s“ 1 for 
CH3O2 to 2 x 10 ~ 26 cm 3 molec“^s” 3 for NO 2 . The formation of H2O- 
SO2 complexes is undoubtedly the basis for the mechanism by which 
H 2 0 enhances the rate of the N 0 2 reaction by many orders of 
magnitude. It appears likely that such complexes may react with 
other oxidants at rates much higher than those with SO2 and that 
one or more such reactions may make a significant contribution to 
the observed rate of S 0 2 oxidation in the atmosphere. 

The dependence of this H2O catalyzed mechanism on oxidant 
concentrations would lead to the observed seasonal dependence of 
SO2 oxidation rates. (1) 

In summary, an energy gradient program was used to locate 
stable structures for H 2 0 • S 0 2 complexes and SCF calculations 
were carried out to determine the binding energies of complexes 
with multiple water molecules. A 4 - 31 G basis set was used for 
most potential energy searches. More accurate basis sets including 
a generally contracted basis set with d orbitals on the sulfur 
atom were used for geometry and binding energy verification. For 
one water complexes five different stable geometries were located 
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with adjusted binding energies betweem 3 and 8 Kcal mol -1 suggesting 
a binding shell for H 2 0 around S0 2 . Very little charge transfer 
between S0 2 and H 2 0 was present. Addition of more than one H 2 0 
was found to be energetically favorable and approximately additive 
up to four added waters. An alternative mechanism for the trospheric, 
gas phase production of acid rain is suggested. A process for 
the formation of an S0 2 containing water droplet is advanced. 
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Table 1. 


SCF total energies (Hartrees) , binding energies (Kcal mol -1 ), in 
parentheses, and total charge on SO 2 of SO 2 (H 2 0) n complexes. 


Geometry Optimized Geometry Fixed 


h 2 0 

-75.9086 

-76.0110 

-76.0110 

-76 . 012 1 

- 

so 2 

-546.3733 

-546.9719 

-547.1540 

-547 . 1755 

- 

I 

-622.2894(4.7) 

- 

- 

-623 . 1930(3.4) 

-0.003 

II 

-622.2943(7.8) 

- 

- 

- 

- 

III 

-622 .2894 (4.7) 

- 

- 

- 

- 

IV 

-622.2902(5.8) 

- 

- 

- 

- 

V 

-622.2993(10.9) 

-622.9974(9.1) 

-623.1770(7.5) 

-623 . 1998 (7.6) 

-0.035 

VI 

-698.2120(13.5) 

- 

- 

- 

- 

VII 

-698.2236(20.8) 

- 

- 

-699 .2212(13. 5) 

-0.067 

VIII 

- 

- 

- 

-775.2437(20.0) 

- 0.061 

IX 

- 

- 

- 

-851.2560(20. 1) 

-0 . 056 

X 

- 

- 

- 

-927 . 2703 (21.5) 

- 

Basic 

Set A 

B 

C 

D 

D 


A: 4-31G 

B: Double Zeta S(12s/8p)/<62211/621>, 0 (9s/5p)/<6111/41>, H(4s)/<31> 

C: Double Zeta + Sulfur d 

D: S ( 12s/8p/ld) /<4s/3p/ld> , 0 (9s/5p) /<3s/2p> , H(4s)/<2s> 


Table 


H 2 0 

S0 2 

I 

II 

III 

IV 
V 

VI 

VII 


2. Bond distances, bond angles, and dipole moments for fragments and 
complexes from SCF geometry optimizations 


Basis 

set 

R(SO) 


L(OSO) 

R (OH) 


L(HOH) 

LA 

LB 

RA Dipole 
moment 

A 

— 


— 

0.951 


111.2 

-- 


2 . 49 

B 

— 


— 

0 .952 


112 . 5 

— 

-- 

2 . 53 

A 

1.530 


114 . 2 

— 


— 

— 

-- 

3 .30 

B 

1.535 


113 . 0 

— 


— 

— 

— 

3 .44 

C 

1.428 


118 . 1 

— 


— 

— 

— 

2 . 67 

A 

1.531 


113 . 1 

0.951 


111.0 

— 

3 . 27 
2.86 

— 

A 

1.528 / 

1.534 

112.9 

0.953 / 

0.949 

112 . 1 

— 

2 . 54 

— — 

A 

1.525 / 

1.531 

114 . 3 

0.953 / 

0.950 

110.6 

— 

2 . 02 

6. 10 

A 

1.531 


112 . 8 

0.951 


110.3 

90.4 

2 . 37 

6.35 

A 

1.533 


113.6 

0.953 


112 . 5 

124.4 

89.4 

2 .48 

2 .95 

B 

1.541 


112 . 6 

0.953 


113 . 2 

125.9 

87.6 

2.49 

2.92 

C 

1.427 


115.7 

0.9 53 


113 . 8 

129.0 

2 . 52 
2 . 85 

2 .51 

A 

1.534 


112 . 5 

0.951 / 

0.949 

112 . 3 

90.3 

2 . 54 


A 

1.541 


113.0 

0.953 


112 . 5 

124.7 

2 .48 

0.07 

All distances 

are in 

Angstroms . 

Dipole 

moments 

are in 

Debyes . 

A second 


set of bond lengths is given for asymmetric complexes. 




FIGURE 1. CONFORMATIONS OF S0 2 H 2 0 COMPLEXES 
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FIGURE 2. CONFORMATIONS OF MULTIPLE WATER 
COMPLEXES OF SO*. 


